o 



Equilibrium relationships for non-equilibrium 
chemical dependencies 

G.S. Yablonsky^'b'* D. Constales^'i G.B. Marin ^ 

^ Parks College, Department of Chemistry, Saint Louis University, 34-50 Lindell 
Blvd, Saint Louis, MO 63103, USA. 

^The Langmuir Research Institute, 106 Crimson Oaks Court, Lake Saint Louis, 
CN M063367, USA 

^ ^Department of Mathematical Analysis, Chent University, Galglaan 2, B-9000 

1"^ Cent, Belgium. 

^^ '^Laboratory for Chemical Technology, Ghent University, Krijgslaan 281 (S5), 

B-9000 Gent, Belgium. 

1-^ Abstract 

G In contrast to common opinion, it is shown that equihbrium constants determine 
the time-dependent behavior of particular ratios of concentrations for any system of 

'""' reversible first-order reactions. Indeed, some special ratios actually coincide with the 

J^ equilibrium constant at any moment in time. This is established for batch reactors, 

■^ and similar relations hold for steady-state plug- flow reactors, replacing astronomic 

^^ time by residence time. Such relationships can be termed time invariants of chemical 

. kinetics. 
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1 Introduction 



In presenting the foundations of physical chemistry, the basic difference be- 
tween equilibrium chemical thermodynamics and chemical kinetics is always 
stressed. 
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A typical problem of equilibrium chemical thermodynamics (ECT) is calcu- 
lating the composition of a chemical mixture that reacts in a closed system 
for an infinitely long time. ECT does not consider time. 

In opposite, chemical kinetics is the science of the evolution of chemical com- 
position in time. 

Some selected results of theoretical chemical kinetics are obtained from ther- 
modynamic principles, especially the principle of detailed equilibrium: 

(a) the uniqueness and stability of the equilibrium in any closed system, see 



(Zeldovich, 1938) and the analysis in (Yablonsky et al. , 1991); 



(b) the absence of damped oscillations near the point of detailed equilibrium. 



see (Wei and Prater, 1962) and the analysis in (Yablonsky et al. , 1991); 



some limitations on the kinetic relaxation from the given initial conditions, 
e.g., based on the known set of equilibrium constants that determine the 
equilibrium composition, one can find a forbidden domain of compositions 



that is impossible to reach from the given boundary conditions. See ( Gorban 



et al., 1982 Gorban, 1984 Gorban, Kaganovich et al. 2006) 



However, the present dogma of physical chemistry holds that it is impossible 
to present an expression for any non-steady state chemical system based on its 
description under equilibrium conditions, except for some relations describing 
the behavior in the linear vicinity of equilibrium. 

The goal of this short note is to announce that, contrary to this 'dogma', we 
have obtained relations of equilibrium type for some non-steady state chemical 
systems. This has been achieved for all linear cases (with a general proof) and 
some non-linear ones. In a forthcoming extended paper these relationships will 
be explained in more detail. 

Based on our results, equilibrium thermodynamic relationships can be consid- 
ered not only as a description of the final point of temporal evolution, but as 
inherent characteristics of the dynamic picture. 

The following are examples to illustrate our statement; in all of them we 
analyze traditional models of chemical kinetics based on the mass-action law. 
The processes described by these models occur in a closed non-steady-state 
chemical system with perfect mixing (batch reactor) or in an open steady-state 
chemical system with no radial gradient (plug-flow reactor), whose description 
is identical to the batch reactor, replacing astronomic time by residence time. 

In this short communication, we shall analyze a combination of data of differ- 
ent thought experiments, e.g. 



a chemical reactor is primed with a substance A only. 



• it is primed with substance B only. 

In both cases, we shall monitor both concentrations A and B and pay special 
attention to the dependencies of "-B produced from A" and "A produced 
from B" . We shall use the notation ^^(t) for the temporal concentration 
dependence of substance A, given the initial condition {A, B, . . .) = (1,0,...), 
i.e., only A occurs, with normalized concentration 1. Similarly, BA{t) is the 
concentration of substance B for the same initial condition; Asit) that of 
A when at t = 0, {A, B, . . .) = (0, 1, . . .). We found interesting relationships 
between these dependencies. 



2 Linear cases 



2.1 A single, first-order reversible reaction A^ B 



By elementary mathematical techniques, for t > 0, 
_ K + kt exp(-(A;+ + kY)t) 

BA{t) 

^-(^) kt + K 

kt + K exp(-(A:+ + ki)t) 

""^^'^^ WTk ■ 



Comparing B^ and Ab, it is clear that they are constantly in a fixed propor- 
tion: 

S.4(i) kt 



which is the thermodynamic constant of equilibrium. It is remarkable that this 
ratio holds at any moment of time t > 0, and not merely in the limit t — )■ +oo. 





K-\ + n,\ 


kt{l - 


- exp{-{kt + k^)t)) 




rC^ -\- rC^ 


A;r(l- 


-exp{-{kt + ki)t)) 



2.2 Two consecutive first-order reactions, the first being reversible: A ^ B ^ 

k- 

C 



Using, e.g., Laplace domain techniques, the following analytical solutions are 
obtained: if we define the expressions 



k+ + k^ + kt ± ^{kt + K + ktY - Aktkt 

-^1,2 — ;^ (2) 



and verify that 

Ai > A;+ > A2 > 0, Ai > A;+ > A2 > 0, (3) 

we can write 

AA{t) = j^rvT ^ (Xi{k+ - A2) exp(-A2t) + A2(Ai - k^) exp(-Ait)) , 

A;2 (Ai - A2) ^ ^ 

k+ 
BaH) = , ' , (exp(-A2t) - exp(-Ait)) , 
M — A2 

CA{t) = 1 - T ^ (Ai exp(-A2t) - A2 exp(-Ait)) , 

Ai — A2 

Anit) = , ' , (exp(-A2t) - exp(-Ait)) , 
Ai — A2 

BBit) = y+TT ^ (A2(Ai - kt) exp(-A2t) + \iik+ - A2) exp(-Ait)) , 

A;2 (Ai - A2) ^ ^ 

Csit) = 1 - Y^-^ ((Ai - kt) exp(-A2t) + {k^ - A2) exp(-Ait)) . 
Ai — A2 ^ ^ 

Again, it turns out that Ba and Ab are always in fixed proportion: 
BA(t) kt 

and again this is the constant of equilibrium of the A ^ B reaction. 



h^ h'^ h^ 

re, K, fc^ 

2.3 The cycle of three reversible first-order reactions A ^ B ^ C ^ A 



fcj fej fcg 



In the Laplace domain, defining the symbols 



(T\ — rb-| ~p nJi ~|~ /Co ~|~ /Co ~p /Co ~p /Co 

02 — /Ci /Co ~p /Co /Co ~p /Co /C-| ~p /Ci /Co ~p /Co /Co "n /Co /Ci ~p 
/Ci /Co '] /Co /C-| ~|~ /Co /Co 

A(s) = s(s^ + ais + 0-2), 
the transformed concentrations A^, Ba and A^ are given by 



CAa{s) 
CBa{s) 
CAb{s) 



S \ t/Ci ~r /Cq ~r /Co ~r /Co jS \ tfCi /Co ~|~ /Ci /Co ~r /Co /Co 

W) 

/Ci 5 ~n I /Ci /Co ~|~ /Ci /Co ~|~ /Co /C"? ) 
/Ci iS ~|~ I /Ci /Co ~|~ /Ci /Co n" /Co /Co ) 



and the ratio of the latter two, by 



CBa{s) kfs + {kfk^ + kfk^ + fc^A;^) 



CAb{s) 


/Ci S ~r \/^i '^'X \ /^i /^2 "1" 2 3 / 
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where the Onsager relationship kfk2k^ = fcf fc^fc^ was used in the final step. 
Hence for all s, CBa{s) and CAb{s) are in fixed proportion given by the 
equilibrium constant (kf/ki). Since the inverse Laplace transform is linear, 
the same proportion holds in the time domain: 



BAJt) ^ kt 



K. 



eq,l' 



(5) 



Similarly, 



rCn 



Beit) k^ 



K. 



eq,2' 



Ac{t) _ kj 
CA{t) k-^ 



K. 



eq,3' 



(6) 



As an example, we show in Fig. [I] the time dependence of Ba/Aa-, Bb/Ab 
and the time-invariant ratio Ba/Ab-, for the case of isomerization of butenes 



in dWei and Praterf |1962D, eq. (129): 



4.623 3.724 3.371 

czs-2-butene :^ 1-butene :^ trans-z-hntene :^ czs-2-butene. 

10.344 1.000 5.616 



This behavior is typical of the examples given here: if the system is started 
from the initial values {A, B,C) = (1, 0, 0), the ratio B^/Aa grows at first from 
0, eventually reaching the equilibrium value (note that for these parameters 
Ba/Aa slightly overshoots the limit). Similarly, when started from (0,1,0), 
the corresponding Bb/Ab initially decreases from +oo and eventually reaches 
same limit, viz the equilibrium ratio. But surprisingly the combination Ba/Ab 
is constantly equal to the equilibrium value, for all times t > 0. 



£■+ h'^ h^ h'^ 

rvj-i ti^cy A>o "'/I 

2.4 The cycle of four reversible first-order reactions A^B^C^D^A 



Although the expressions become more involved, it is straightforward to verify 
that in this case the fixed equilibrium proportions also hold, given the Onsager 
relationship kik2k^kf = kik2k^k^: 

BaH) _kl_ 

Asit) ~ kl -""eq,!, UJ 



and similarly for Cb/Bq, Dc/Cd and Ad/ Da- Furthermore 



CAJt) ^ k+kt 

Acit) k^k^ 



^eq,l^eq,2> 



with a similar relationship for Db/Bd. 

2.5 Sketch of proof for general systems of first- order reactions 



We use the terminology and results of (Roelant et al. , 2010). Let A and B 
be substances such that a path from either to the other exists. In view of 
Onsager relations, then there must exist a reversible path from A to B. The 
denominators of Ab and Ba in the Laplace domain are the same, but their 
numerators differ. For Ab-, the numerator consists of contributions from all 
forests where A is the root of a tree and B is in that tree. Let Z, X, Y denote 



other nodes as in the denominator of the left-hand term in ([9]): 
Z ^ A^ X ^ B ^Y kA^xkx^B 



Z ^ A^X ^B ^Y k 



K 



A^. 



>B 



(9) 



X^AI^B^X 



Reversing the X -^ A and B -^ X arrows amounts to muhiplying the forest's 
term by (kA^xkx^s) / (kx^Aks^x) ■ In view of the Onsager relations, this 
value does not depend on X, but only on A and B; in fact, it is the equilibrium 
constant K^^ ^^^b- If several X or Y or Z occur, the reasoning is the same. 
Consequently, every term in the numerator of Ab is in that proportion to the 
corresponding term in the numerator of Ba, and the fixed proportion for all 
s translates directly to the time domain: 



ABit) 



K, 



¥B- 



(10) 



3 Nonlinear cases 



3.1 Nonlinear reversible reaction (forward second order, backward first order) 

2AtB 

k- 



In accordance to the mass conservation law, the balance is 



A{t) + 2B{t) = l 

In order for the A and B trajectories to reach the same equilibrium, we choose 
to start them from (1, 0) and (0, 1/2) respectively. The nonlinear differential 
equation is 



dA{t) 
dt 



-2ktA^{t) + k:[{l-A{t)) 



which can be solved analytically as 



AA{t) 



. 8-^ + 1 + tanh -tkf 
\ fci \ 2 






+ 1 



^.| + l+(4| + l)tanh it.- 






BA{t) 



2-^tanh -tkt 
K \2 ' 






.8^ + 1+^4^ + ljtanh -t.- 



M 



Asit) 



2 tanh | -tk'l 






+ 1 



A;+ /I 

. op + l + tanh -tA;+ 






The remarkable proportion in this case differs shghtly from the hnear exam- 
ples: 



K^ 



Ba _ kt _ 
AaAb ~ ki - ""^^'l 

we see that the denominator involves the A concentrations of both trajectories, 
A A and Ab- Only this can ensure a ratio that equals the equilibrium constant 
at every time t > 0. 



3.2 Nonlinear reversible reaction (forward and backward second order) 2 A ^ 2B 

k- 



The mass conservation law is 



A{t) + B{t) = l, 

which offers no difficulties for the initial values, (1, 0) and (0, 1). The differen- 
tial equation 



dA{t) 
dt 



2ktA\t) + 2k^{l-A{t)y 



can be solved analytically as 



AA{t) 



1 + 



f4 



tanh 2tJk{kl 



BAit) 






tanh 2tJk{k'^ 



1 + 



kt 



\K 



'- tanh (2t J k{kl 



Anit) 



Bsit) 






tanh 2tJk{k^ 



1 + 



ki 

\ik 



^ tanh \2t\lktK) 



ki_ 

\kt 



tanh 2tJk{kl 



Again ehminating time, the similar proportion for this case is 



BaBb ^ kt 

AaAb k^ 



K. 



eq,l' 



where both numerator and denominator have undergone a duphcation in A 
and B trajectories, still producing the equilibrium constant at all times t > 0. 



4 Conclusions 



We are going to describe these results in detail in a full-length paper. Then 
the similar approach will be applied to different systems, i.e., CSTR reactors, 
to reaction-diffusion TAP systems, etc. 

Presenting simply the result of this paper, it is a surprising relationship be- 
tween A from B and B from A. An "ABBA rule" for short, which reminds 
the authors of their youth. 

"They speak strangely hut I understand" — ABBA, "Eagle" (1978). 
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Caption to figure 



B_A/A_A 
B_B/A_B 
B A/A B 



(B/A) starting from (A,B)=(0,1) 



0.6 



S B starting from (1 ,0) divided by A starting from (0,1 ) 



0.4 - 



0.2 - 



(B/A) starting from (A,B)=(1 ,0) 



0.05 0.1 



Fig. 1. Time dependence of Ba/Aa, Bb/Ab and the time-invariant ratio Ba/Ab, 



for the case of isomerization of butenes in (Wei and Prater, 1962), eq. (129) 
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